Reconstruction of the bubble nucleating potential 



Luca Amendola^, Carlo Baccigalupi^'^, 
Rostislav Konoplich^, Franco Occhionero^ & Sergei Rubin^ 
^ Osservatorio Astronomico di Roma, Viale del Parco Mellini, 84, 00136 Roma, Raly 
Dipartimento di Fisica, Universitd di Ferrara, Via Savonarola 9, 44^00 Ferrara, Raly 
'^Moscow Engineering Physics Institute, Kashirskoe sh, 31, 115409 Moscow, Russia 



Abstract 



We calculate analytically the bubble nucleation rate in a model of first 
order inflation which is able to produce large scale structure. The compu- 
tation includes the first-order departure from the thin-wall limit, the explicit 
derivation of the pre-exponential factor, and the gravitational correction. The 
resulting bubble spectrum is then compared with constraints from the large 
scale structure and the microwave background. We show that there are models 
which pass all the constraints and produce bubble-like perturbations of inter- 
esting size. Furthermore, we show that it is in principle possible to reconstruct 
completely the inflationary two-field potential from observations. 

PACS: 98.80 Cq, 98.80 Es 
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I. INTRODUCTION. 



One of the most interesting ideas introduced in inflationary cosmology in recent years is 
the possibility of performing a phase transition during inflation. In such scenarios, two fields 
act on stage: one, say uj, slow rolls, driving enough inflation to solve the standard problems; 
the second field, say tunnels from a false vacuum state to an energetically favoured true 
vacuum state, producing bubbles of the new phase embedded in the old one. Both processes 
are governed by a two-field potential U{u!,ip). To avoid the graceful exit problem, the true 
vacuum state has to allow for a period of inflation on its own. We can then speak of a true 
vacuum channel over which the bubbles slow roll until inflation ends, and reheating takes 
over. Depending on the potential, three classes of first-order inflation models have been 
proposed so far. The first is the classical extended inflation the bubbles are produced 

in a copious quantity, so that they eventually fill the space and complete the transition. 
To avoid too large distortions on the CMB, this scenario must produce very small bubbles 
0J^, so that they are rapidly thermalized after inflation. No trace of the bubbles is left 
in our Universe, and from this point of view such scenarios do not lead to new predictions 
over inflation without bubble production. The second class is the < 1 inflation 
here the transition is never completed, so that each bubble resembles an open Universe to 
inside observers. Therefore, if the bubbles inflate for less than the canonical A^^ = 60 or so 
e-foldings, they will approach an Q < 1 Universe. Here the effect of the nucleation process 
is observable, although it is indistinguishable from a slow roll inflation just shorter than 
Nt e-foldings and with no first-order phase transition at all. Finally, in a third class 
of models has been proposed, following an early suggestion of La [|l^]. In such models, the 
phase transition is completed before the end of inflation. Then, it has been shown that the 
primordial bubbles can be large enough to drive structure formation, and still be below the 



CMB level of detection ||TTHl^. In such a scenario, the present large scale structure is a 



direct outcome of the first order transition, which is therefore observable and testable. 



Several deep redshift surveys detected large voids in the galaxy distribution ||T5|-|T7 
although it is still not clear if they are really empty of matter or just lack luminous galaxies 
||18||. Standard models of galaxy formation can barely accounts for these structure, and do so 



only at the price of adjusting the parameters to get very large scale power (see e.g. Ref. ||T9|). 
Therefore, just as we associate matter clumps to primordial fluctuations, it appears worth 
trying to associate the present voids to primordial bubble-like fluctuations, produced during 
a first-order phase transition. Within different contexts, the idea of the voids as separate 



dynamical entities has been investigated several times in earlier literature pO |. 

A crucial aspect of bubble inflation is the calculation of the bubble spectrum nsiL), 
defined as the number of bubbles per horizon with comoving size larger than L. In M we 



calculated in a specific model, built on fourth order gravity which we found to 
possess the requested features. We found that ub^L) can be approximated by a power law, 

nB = {L^/LY, (1) 

and that Lm can be as large as the observed voids in the Universe. 

The central quantity needed to evaluate is the nucleation rate in the semi classical 
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r = A^^exp(-5) (2) 
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where is a constant with a dimension of mass, and B is the Euchdean least action 
minus the action for the external deSitter spacetime solution. The calculation of F in |Q 
(and in most other papers on the topic) was based on the thin-wall limit, on neglecting the 
gravitational correction, and on a dimensional argument for evaluating Ai. In this paper we 
remove, at least partially, all three approximations. This calculation allows us to reconstruct 
completely the inflationary two-field potential from the determination of four observable 
quantities: the slope of the bubble spectrum, its amplitude, the density contrast inside the 
bubbles, and the amplitude of the ordinary slow-rolling fluctuations. We remark that only in 
a model, like ours, in which the bubbles are directly observable, it is possible to reconstruct 
the tunneling sector of the primordial potential. 

The scheme of the paper is as follows. In Sec. 2 we present the class of potentials we 
are going to investigate. In Sec. 3 we give the basic formulas. In Sec. 4,5, and 6 we go 
through the detailed calculation of F, taking into account deviations from the thin wall limit, 
inserting the gravitational correction and calculating explicitely the factor Ai. In Sec. 7 we 
introduce our model of first-order inflation in a fourth-order gravity theory and calculate 
the time-dependent nucleation rate. In Sec. 8 we write down the bubble spectrum and in 
Sec. 9 finally we compare it with constraints from CMB and large scale structure. In the 
last section we draw our conclusions. 



II. THE MODEL 

We consider the scalar field theory described by the action (hereinafter, h = c = 1) 

S = /A^{-^ + lvV#'--CW} (3) 

where g is the metric determinant and TZ the curvature scalar. The potential is a generic 
quartic function with non-degenerate minima which allows for tunneling. We can write it 
very generally in the form 

U{ij)=A + Vi{^) + V2W (4) 

where A is a cosmo logical constant, Vi is a quartic with two equal-energy minima and V2 is 
a symmetry-breaking potential which brings the energy of one minimum, the false vacuum 
(subscript F), to a value larger than the other, the true vacuum (subscript T). In the two 
minima, the Einstein equation reduces simply to 

= 87VGU/3 (5) 

We will denote with Ht, Ut the Hubble constant and the potential energy of the true vacuum, 
and with Hp^ Up the same quantities for the false vacuum. We wish to calculate the tunneling 
rate (0) where B = Spijp) — Sp^ipp), and Sp^ip) is the Euclidean least action, i.e. the Action 
of the "bounce" solution. We perform the computation in the thin-wall limit (actually we go 
to a post-thin- wall limit calculation), according to which the 0(4) bubbles nucleated have 
4-radius R ^ A, where A is the wall thickness. Further, we include the gravitational term 
in the action: as we will show, this term is important when the parameter g = RHt (not to 
be confused with the metric determinant) is much larger than unity. This limit amounts in 
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fact to a bubble approaching the space curvature radius 1/H. It is convenient to write the 
potential @ in the form 

Then, the true vacuum state is ip = —ipo, and the false vacuum is ip = ipQ. The potential (|) 
is therefore defined by four physical parameters: g, R, A, -00- 



III. THIN WALL LIMIT 

The Euclidean action of the scalar theory is 

7^ 1 



Se = J d^x^^^ 



167rG T'^^ 



(7) 



In the Euclidean metric for a 0(4) space ds^ = dr"^ + a'^{r)dQ'^ one has 7^ = —6{aa" + a'^ 
l)/a^ and 



Se = 27T^ / dr 



-^{a^a" + aa'^ - a) + a^^^ + U) 
ottG 2 



where the prime denotes derivation with respect to the 4-radius r. The Euclidean Klein- 
Gordon equation for ip is 



ip" + 3-ij' = dU/dip 
a 

and the Euclidean Friedmann equation is 



j2 



1 + - u) 



(9) 



(10) 



In the zero-gravity limit, G — 0, the latter equation gives a' = const, so that (^ reduces to 

(11) 



iP" + -t/j' = dU/dip 
r 

In the thin- wall limit, in which A, one can assume that the second term in (pUD can be 
neglected, and that dU/dip = 2%p{'il)'^ /ip^ — 1)/A^. The solution which interpolates between 
false and true vacuum is then 



ip^^'^ = ipo tanh 



r - R^ 
A 



(12) 



where Rw is an integration constant that will be determined later. To integrate the action 
over the bounce solution, we consider that outside the bubble, i.e. in the false vacuum, 
ip = ipo, so that Bext = SEii'o) — SEii'o) = 0. On the wall, at distance R^, we have 



B 



wall 



2'k'^R^Si 



(13) 
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where 



^1 



-^0 



(14) 



Finally, inside the bubble, ^p = —Tpo, and since (from (p!OD, and neglecting ■?/'') da 
dr [1 — a?H'^^^'^ we have 



B 



int 



2G 



ada 



[l-a^H^^fl^-il-a'HlY/^ 



(15) 



The general expression is therefore 



B{R^) = 2n'RlS, + ^ {h^' [(1 - RlH^f' - 1 



[l-RlHlf'-l]] (16) 



Let us note that R = 3Si/e, where e = Up — Ut = 4:iIjq/RA. Then we see that B{Rw) is 



minimized by 



R^ = R(l+g^ + AnGRSi + Un^G^R^Sl 



-1/2 



(17) 



Then, for G ^ 0, which implies g = RHt 0, one has R^ = R- Notice that the parameter 
RwHt which appears in ([TBI) equals g, since = g^/R^ for the potential @. This shows 
explicitely the role played by the constants R and g. Finally, we obtain the usual (zero- 
gravity, thin-wall) result ||2l| 

(18) 



IV. THE GRAVITATIONAL CORRECTION 



Now to the general case, G ^ 0. To the purpose of this paper, we simplify our problem 
by putting [Hp — H^)/H^ <C 1, i.e. assuming that the vacuum energy difference is much 
smaller than the true vacuum energy. Our results will be consistent with this approximation. 
This is equivalent to neglecting the last two terms in parentheses in Eq. ([T7|) . Then we have 
that B{R^) is minimized by a bubble radius R^ = -R/(l + g'^Y^'^- The bounce action is 
B = B^f^g) where Bq is the no-gravity action (|18]) and where 



f{g) = 4(1 + /)-3/2 {l + g-'[2 + 3g' - 2(1 + gY']} (19) 
For g ^ 0, f{g) 1, as expected. To the lowest non-trivial order in g, 

B = Boil-g^) (20) 



Gravitational effects, then, increase the nucleation rate F pT . 
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V. POST-THIN- WALL CORRECTION 



Let us come back to the Euclidean spherically symmetric equation of motion in absence 
of gravity, Eq. (|TT]). While the trivial true- vacuum solution ■?/'_ = —ipQ holds also when the 
symmetry-breaking term in the potential is considered, the false vacuum solution is slightly 
displaced: 



/ A A2 ^ 



(21) 



It is useful now to expand the equation of motion around r = R. Introducing the variable 
z = {r — R)/A we obtain 



(22) 



where now the prime denotes derivation with respect to z. We search solutions of Eq. (^ 
to the first order in A/R: 



^ = +^(i)(^/2i?) + ... 



(23) 



We know already that to the zero-th order ^^'^^ = ipotanhz. Subtracting the zero-th solution 
from Eq. ( ^2]) we get, for i > 1 



-4 



6 



cosh z. 



(24) 



where fi{z) is defined order by order from Eq. (^) and depends only on z and not on 7/'*^*^ 
At first order, then, the solution is 



iIj = ^lJo{t&nhz + ( 



cosh z 



— tanh z — 1 



2(^) [tanh2;+l] 



(25) 



It can be noticed, from Fig. 1, that the post-thin- wall correction moves the large-z value of 
i/j from ipo to ?/'+. In Eq.(p5D we included the part of the correction at second order which 
also contributes to the leading classical correction to the action after integration over z. 
The Euclidean action, finally, is calculated as 



B = S{^) - Siipo) = Bo{l - 9A/2R) 



(26) 



The term 9A/2R is therefore the post-thin- wall correction term, which again increases the 
tunneling rate. 



VI. THE PRE-EXPONENTIAL FACTOR 



In semiclassical approximation in four dimensions the rate of false vacuum decay is given 



by Coleman's formula pT 



6 



r 



2-k) 



det^ 



-d^ + U" (^)) 



det (-92 + U" 



-1/2 



(27) 



where the functional determinant is computed with the four zero eigenvalues omitted, if) is 
the classical solution of Eq. ([2^), ipj^ is the trivial false vacuum solution. It was shown in 
that in the thin wall limit the functional determinant can be expanded in power series 
of 2i?/A as 



In 



det^^^ D 
det(^) Do 



Ci ln(2i?/A) + C3 {2R/d.f + C2 {2R/Af + .. 



(28) 



The parts of this expression, which behave as powers of 2R/A have no universal significance, 
for they are adjusted or even completely concealed by ultraviolet renormalisations and we 
will include them in overall multiplicative factor A4'. The logarithmic part, on the other 
hand, is not expected to be affected by short distance structure. The universal infrared 
logarithm ln(2i?/A) is associated with the geometrical features of the bubble and does not 
depend on the specific form of the double well potential. Since the logarithmic term appears 
due to the infrared region we have to consider the low-lying levels of the operator of the 
small fluctuations around the bounce solution 



D 



+ U" iij) . 



(29) 



This operator is rotationally invariant and thus its eigenfunctions are, in their angular de- 
pendence, 4-dimensional scalar spherical harmonics. In terms of radial eigenfunctions the 
eigenfunction equation becomes 



dr 



fni{r) = Kifni{r) 



(30) 



Expanding the centrifugal potential in power series around r = i? (it is justified when we 
deal with the low lying states / -C 2R/A since bound states are localized at the wall of the 
bubble), keeping only the leading term, and replacing r with z, we can rewrite Eq. ( pOD for 

n = as 



d^ 



dz^ 



+ A^U" (^) - 2Eo 



foi = 0, 



(31) 



where 



En = — 
" 2 



A'Ao/ - 4 



4/(/ + 2) + 3' 



(2i?/A)2 

does not depend on the quantum number /. Taking in the form of expansion 



(32) 



Ea 



a + 



+ 



2R/A {2R/A)\ 
we obtain the spectrum of the eigenvalues in the following form 



(33) 
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^ _ 4 
- ^ 



b c 4/(/ + 2) + 3 

^ 2R/A ^ (2i?/A)2 + ^ 4(2/?/A)2 



(34) 



We know that in the spectrum of the eigenvalues there is one negative eigenvalue Aqo corre- 
sponding to uniform expansion of the bubble, which is just vacuum instability. There is also 
the 4- fold degenerate eigenvalue Aio = describing translations of the bubble center in the 
Euclidean space in four directions. Therefore, substituting / = 1 into Eq.(0) we obtain 

a = -1,6 = 0,c= -15/4 (35) 

and find the lowest band of the eigenvalues 

(/-l)(/ + 3) 

= ^2 (36) 

It is important to know that in deriving Eq. (^) for Aq/ we assumed the existence of 
nontrivial classical solution of the equation of motion which break a translational invariance 
(zero modes) and the thin wall limit in order to justify the expansion in powers of 2R/A. 
Therefore this equation does not depend on specific form of the potential U (ip) . 

Calculating the product of Aq; for / -C 2R/A and keeping only the logarithmic term it 
is strightforward to show that q = —10 (see Eq. |28|). Then it follows immediately from 
the general formula (^) that in the thin wall limit, and taking into account the quantum 
corrections, the rate of false vacuum decay is given by 



{2R/A) 

independently of the specific form of the potential U (ip) . 



e-^ (37) 



VII. A FOURTH-ORDER GRAVITY, FIRST-ORDER INFLATION MODEL 



As we mentioned in the Introduction, astrophysically interesting bubbles can be produced 
only in a two field model, in which one has slow-rolling along a field, u, while the second 
field, ip, performs a phase transition in an orthogonal direction. 

We already introduced in a specific model of first-order inflation able to generate large 
bubbles in fourth order gravity. Once a conformal transformation is applied, the model is 
described by the action (we put hereinafter also G = 1 ) 



S 



'—gd^x 



n 

IGvr 



(3/47r)^.^a;''^ + 



1 



-2uj 



where the potential is 



-4aj 



3/lf2 
oZn 



2u\2 



(38) 



(39) 



and where 
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2 ;2 

-m ip . 



(40) 



The slow roll inflation driven by uj takes place at ^ 1, and is over when uj approaches zero. 
At large uj, the potential U is dominated by W{tp), and thus the false vacuum minimum at 
ip ^ ipo, for which Up = e-^'^[(3MV327r)(l - e^^)^ + Viipo)] is unstable with respect to 
tunneling toward the true vacuum ip = (for which Ut ~ 3M^/327r). At small u, U 
is dominated by V^ip), and both the true and the false vacua converge to the global zero- 
energy minimum a.t uj = = 0, where inflation ends and reheating takes place. The slow-roll 
solution in this model for u ^ 1 can be written very conveniently as |^ 



2uj 



(41) 



where A^ is the number of e-foldings to the end of inflation. In the same limit, H = M/2. 
The potential (|39D for 1 takes the form (||) by the substitution: 



^ ^ ^ + Vo/2 
and therefore in the new notations we have: 



R = (3/2 



TT 



,1/2 



A = (87r/3) 
g = MR/2 



1/2 ^0 

MAi/2 



(42) 



(43) 

(44) 
(45) 



Notice that now R depends on the e-folding time as i? ~ A^^. This is where the scale- 
dependence of the bubble spectrum arises from, as we show below. 
The Euclidean action is 



lOTT 2 



(46) 



(neglecting the kinetic energy of uj). To obtain a canonical kinetic term, we rescale the 
coordinates |25] as = e~'^x^ , so that, in the new coordinates 



lOTT 2 



(47) 



where finally Se is canonical. Combining Eq. (^7]), Eq. ( PB| ) and Eq. (^), we obtain 



B 



6 A 



9A 
2RJ 



(48) 



(there is an extra factor of 1/4 because now the vacuum states are at and V'o instead of at 
ib'i/'o)- In terms of A^ this is 



, N2 



41 



(49) 
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where 



N 



33/2 ^pQm^ 

4 M2A 
27tc i)lm^ 

8 M^A 
V 32 j M2AV2 



(50) 
(51) 
(52) 



The two approximations we adopted are vahd for N > N2 (thin- wall), and for N < 
(gravitational correction). It is useful now to introduce a fourth e-folding epoch, Nq, defined 
as the epoch at which the crucial quantity Q = AirT/dH'^ equals unity. This can be seen as 
the epoch at which one bubble for horizon volume for Hubble time (or 4-horizon) is nucleated, 
that is, the bubbles are saturating the false vacuum space. The bulk of the nucleation occurs 
therefore just before A^o^ ans spans only few e-folding times. We can therefore approximate 
Q{N) as 



Q 



where 



AttT 
9H^ 



exp 



Nf 



N2 
No 



No' 
N^ 



In 



[l-{N,/Nom-{No/N,)^] \ 9M4 



(53) 



(54) 



From Eq. (|7|) we can write the factor M as M' B^/'^{A/2R). 

The scalar and gravitational zero-point fluctuations generated during the cj-driven slow- 
rolling are proportional to M |2^. From the current microwave background measurements. 



we obtain M ^ 5 ■ 10 (in Planck units), a value that we will adopt in the numerical 
results below. However, since in our model one should also consider the contribution of 
the bubbles to the microwave background, this constraint is actually only an upper limit 
on M. Once M is fixed, the three e-folding constants A'^i-s fully determine the inflationary 
potential, through the three remaining constants m,\,ipo. A successful sequence of epochs 
is the following: at N^, the gravitational correction ceases to be important; at A^o the crucial 
parameter Q is of order unity, and the bubbles saturate the space; at N2 the thin wall 
approximation breaks down; and finally at A^i the bounce action B equals unity, and thus 
the semi-classical approximation is no longer reliable. For traces of the phase transition to 
be visible today, we must also impose A^o < Nt, being Nt the observable last 60 e-foldings 
of the inflation. Then, the full set of constraints is 



(55) 



VIII. THE BUBBLE SPECTRUM 



We can now calculate explicitely the bubble spectrum in our model, 
bubbles nucleated in the interval dt is 



The number of 



10 



duE 



raV,„exp(-/) , I=\ 1^ dt'T{t')a%t') (£ 



f* dr 




(56) 



where Vin is the horizon volume at = A^^, Vm = Air/SHf^, and where the exponential factor 
accounts for the fraction of space which remains in the false vacuum. To get a manageable 
expression, we first change variable in Eq. (|56|) from the nucleation epoch t to the scale L in 
horizon-crossing at t, by use of the relation dL/dt ^ —HinLh/a valid during slow roll. This 
gives 



dn 



B 



dL 



3LiL^^Q{N)e- 



We approximate Q{N) around N = Nq as Q ^ exp[s(A^o ~ N)], where 

s = 4 — — 



(57) 



(58) 



We make use of the relation between the e-folding time N and the bubble comoving size L 

HL{N) = H.nLh exp{N - Nt) . (59) 

It follows Q = e~'^^^ {Lh/ Ly 1 where AiV = Nt — Nq corresponds to the duration of the 
transition. For / <^ 1, i.e. far from the end of the transition, we obtain 

dUB 



dL 



AL 



-4-s 



and 



where 



Ub = {Lra/LY , p=3 + S 



Lhe 



(3-p)AAf/p 



(3/p) 



i/p 



(60) 



(61) 



(62) 



In the thin- wall, zero-gravity limit, s = AN^/Nf [y]. When / approaches unity, the nucleation 
process reaches a peak, and then declines rapidly, due to the fast decrease in the false vacuum 
space available. In we showed that the peak occurs, as expected, just after Nq, which we 
therefore consider the end of the nucleation. 

As we discussed in Ref. 0, Eq. (^) differs from the analogous quantity in extended 
inflation |T| ||^ |TD| because here we have AA^ instead of A^t ~ 60. We show below, by a 
suitable choice of AA^, i.e. of the duration of the transition, that we may have a bubble 
spectrum hb which generates astrophysically large bubbles and escapes the CMB bounds, 
contrary to what occurs in the other models of first-order inflation. 



IX. RESULTS AND DISCUSSION. 

We wish now to compare the bubble spectrum with the constraints from large-scale 
structure and the cosmic microwave background. It is not the aim of this work to determine 
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the best model parameters, also because the statistics on observed voids is still at a very 
schematic level. We will focus instead on showing that there is a region of parametric space 
which gives astrophysically large bubbles, that can contribute significantly, if not exclusively, 
to structure formation, while passing the CMB constraints. 

When a bubble nucleates, its density contrast can be estimated (in the thin- wall limit), 

as 

6 ^ \6p/p\ = = [(iv/iv,)2 + (63) 

Up 

where 

r^l = (64) 

Since the bubbles cross out the horizon soon after their nucleation, the same density contrast 
will be found at reenter. The bubble density contrast increases to unity for N -C A^4, i.e. 
toward the end of inflation. One can also write iV4 1.5MN^ /N^. For instance, taking 
acceptable values as = 90, A^i = 20 and M = 5 ■ 10~^, one obtains 6^6- 10~^ at 
iV iVo ~ 50. Taking N-^ oo one would obtain completely void bubbles, 6 1. Here we 



will consider instead only small values of 6 [|T3|. From Eq. (|59|) and ^ A^4, we have the 
density spectrum at reenter 



6{L) 



N4_ 



'^-^log(L/L.) 



(65) 



Since 5 <^ 1, we do not consider (contrary to Ref. |^,rT,n]) the overcomoving growth that 
takes place only if and when, in its late history, the density becomes non-linear. 

Let us discuss first the CMB constraints. A bubble of radius L at decoupling produces a 
Sachs- Wolfe distortion on the microwave temperature of AT/T ~ S{L)L'^ /L"^, if denotes 
the horizon scale at decoupling. In reality, bubbles which reenter before decoupling have time 
to deepen; for the range of scales we are interested in, however, this is a minor effect and 
we neglect it. In a pixel corresponding to a size of Lp > L at decoupling, a further factor of 
L'^/Lp smears the signal There are two main CMB constraints arising from observational 
upper bounds to such Sachs- Wolfe effect. Full-sky, low-resolution surveys like COBE can 
detect rare big bubbles as hot spots. On the other hand, a large number of small bubbles 
can be detected as Poissonian fluctuations in high resolution, small coverage experiments 
with antenna beam below l". Here we restrict our attention to the first kind of observational 
effect, as the second one depends on the lesser known physics of the small bubbles nucleated 
near the end of the first order transition. At any rate, the small bubble constraint is expected 



to be less severe than the hot-spot test for bubble spectrum slopes p < 10 |[TT[]. In ||T3| , p!^ we 
also analyzed the constraints on bubble models from Doppler effects on the last scattering 
surface, and found them to be generally smaller than the Sachs- Wolfe ones. 

Assuming a power-law spectrum like (^) the constraint can be put in form of restrictions 
on the two parameters p and for large p [|1^. Let denote the smallest bubble at 
decoupling that can give an observable 3a signal (AT/T ^ 5 ■ 10~^) in a COBE pixel. We 
simply require that there are fewer than one bubble larger than intersecting the last 
scattering surface. Then, the hot-spot constraint amounts to |^J^ 
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1/p 



(66) 



L„ is calculated as 



AT/T 



3 L^L"^ 



(67) 



We will take Ld = 190h^^ Mpc, and Lp = 300h^^ Mpc. The latter value corresponds to a 
beam angular opening of 3'^, which is roughly the COBE beam opening. 

Next we impose to our bubble spectrum the qualitative requirement that it be able to 
produce a significant large scale structure. We realize this condition in a very simplified way: 
we just find the parametric region for which bubbles larger than 3h~^ Mpc fill the space for 
more than 50% a constraint which can be compared, for instance, with the observed 
voids in the SSRS2 database |]T^. Notwithstanding its simplicity, this minimal condition 
puts a strong constraint on the parameter space. If we had indications on the real matter 
content of large voids, we could put a further direct constraint on S. For instance, if the 
voids of 20h~^ Mpc (which are therefore nucleated 5-6 e-foldings after Nt) that are currently 
detected on the large scale surveys (see, e.g., fl^ ) are completely empty of matter, they 
should have 6 > 10~^ at decoupling, in order for them to be cleared of matter by linear 
growth at the present. If some matter, maybe in the form of Lya clouds [jl8|, is found inside 
the voids, then one should impose on the contrary S < 10~^. Here, to fix the ideas, we put 
6 = 10~^ on scales of ~ 30h~^ Mpc. In any case, the results are not very sensitive to 6. 

The main results are contained in Fig. 2. On the plane {p,Lfn), we display as a shaded 
area the parametric region of cosmological interest, i.e. the models which satisfy the CMB 
constraint (0), and fill the space as requested. The CMB constraint is plotted for N3 = 90 
and 6 = 10~^, which, along with M = 5 ■ 10~®, implies A^i ~ 17. These values satisfy 
all the conditions for a successful inflation. Then we plot two curves Lm{p) given in Eq. 
(|6T|), for A^'o = 53,55. It can be seen that these curves cross the acceptable region; our 
model is therefore capable of producing bubble spectra which pass the CMB tests and have 
interesting large scale features. Let us consider one pair of parameters for reference: p = 9, 
Lm = lOO/i"^ Mpc, which lies on the curve A^^q = 54, and is inside the acceptable region. 
Then, from (p5D, and putting 5 = 0.001 and M = 5 ■ 10"^, we have that the full set of 
e-folding parameters is 



No = 54:,Ni = 17,N2 = 16,N3 = 90 



(68) 



which verifies (|55|); we have also N4 = 1.7. Clearly, this is only one of an infinite set of 
parameters that fulfills all the conditions; however, the values are in general close to those 
quoted. We may also note that with such a bubble spectrum, p = 9 and Lm = I00h~^ Mpc, 
one can estimate {Rg/ Lh)^nB{> 20h^^ Mpc) ~ 70 voids larger than 20 h^^ Mpc by radius 
in a spherical survey of depth Rg = 200h^^ Mpc; this value, which of course is extremely 
approximate, is indicative that many, if not all, of the present voids can be explained by 
a first-order phase transition. Finally, the potential parameters, other than M, are (all in 
Planck units) 



A = 4 - 10" 



, m 



3.3- 10~^V^o = 1.8-10" 



(69) 
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Thus, as anticipated in the Introduction, the observations of the two bubble spectrum pa- 
rameters, Lm and p, of the matter content inside bubbles (so to fix 6) , and of the ordinary 
fluctuations driven by u on the microwave background, which give M, are in principle all 
what we need to completely reconstruct the primordial first-order potential. This would 
be impossible in the other models of first-order inflation, in which bubbles do not leave 
observable traces. 



X. CONCLUSIONS 

Large voids are ubiquitous in the Universe. They occupy most of the observable volume, 
and are probably the dominant contribution to the large scale power spectrum [0. The 



standard possibility is that they simply derive from the primordial underdensities which are 
expected in the ordinary Gaussian models of structure formation. However, this scenario 



would be put in jeopardy if larger and larger voids continue to be discovered Moreover, 
if the voids are actually filled with unclustered, or mildly clustered, matter, it would be 
impossible to explain their roughly spherical shape: an underdensity, in fact, becomes more 
and more spherical only as it becomes non-linear pOj. It is therefore worthy to investigate al- 



ternatives. We proposed here, basing on earlier work |]TT| p, a first-order inflationary model 
which produces primordial non-empty bubbles, along with ordinary slow-rolling fluctuations. 
We calculated in detail the nucleation rate, including classical, quantum and gravitational 
corrections, and showed that the model gives a strong contribution to large scale structure, 
while passing the microwave constraints. The determination of four observable quantities 
completely fixes the primordial potential, including the tunneling sector, which, on the con- 
trary, is unobservable in the other models of first-order inflation. 
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FIGURE CAPTION 



• Fig. 1 

The bounce solution, interpolating between false vacuum (FV) and true vacuum (TV). 
The solid line corresponds to A/i? = 0, the dotted line to A/i? = 0.1 and the dashed 
line to A/i? = 0.3. 

• Fig. 2 

Parametric region of the model. The shaded area is the region which passes the con- 
straint from the microwave background isotropy (solid line labelled CMB) and the 
constraint from large scale structure (dashed line labelled LSS). The two dotted hnes 
give L,n = L„i{p) in our model for the two values of A^o shown. The other parameters 
adopted are indicated in the box. 
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